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Abstract 

The electromagnetic excitation of the two-phonon isovector giant dipole res- 



onance in relativistic projectiles incident on heavy targets can be proceed via 
several intermediate 1~ one-phonon giant resonance states. In two step elec- 
ON ' trie dipole transitions the population of + , 1 + , and 2 + two-phonon states 



are possible. We calculate the amplitude distribution of 1~ excitations with 
an RPA formalism, and use it to calculate the electromagnetic excitation of 
two-phonon states in second order perturbation theory and coupled-channels. 
We show that a conspiracy between angular momentum coupling and the 
^ ■ strength of the electromagnetic fields suppresses contributions of 1 + states to 

the total cross sections. 
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The electromagnetic excitation of two-phonon states of giant dipole resonances in heavy- 
ion collisions at relativistic energies has been studied in several experiments fl|-|5| . In LAND 
measurements |T],0J§] these cross sections were extracted from total cross sections by sep- 
arating a contribution coming from excitation of single dipole and quadrupole resonances. 
The observed cross sections are large, on the order of several hundred of millibarn. How- 
ever, the first experimental results for the double giant dipole resonance (DGDR) excitation 
were found to be enhanced by factors of 2-3 for 136 Xe M and 197 Au as compared to 
any theoretical calculations available. Although for the later measurements in 208 Pb || the 
experiment-theory correspondence is much better, theoretical calculations still underpredict 
the DGDR cross section by about 30 %. Another part of the problem with double resonances 
is related to the experimental position and width of these new resonances and how they cor- 
respond to a harmonic picture of nuclear excitation. For more details of experimental studies 
of this problem we refer to reviews ||. 

There are only a few theoretical papers written on this subject. In the first group the 
GDR is treated phenomenologically as a single state and the DGDR is either a sum of these 
two GDR [0 or a sum of two "GDR" phenomenologically mixed with "GDR^" states ||. 
These papers deal mainly with the problem related to the excitation cross section of the 
DGDR. The position and width of the DGDR is a subject for the second group PHT3I where 
the GDR is considered microscopically. In ref. [T(J a fine structure of the GDR was calculated 
by coupling of RPA 1~ states to two-phonon configurations and the DGDR was treated as 
two independent excitations of the single GDR. In other papers eigen wavefunctions of the 
DGDR were used and two-phonon DGDR states were coupled either (a) between themselves 



in ref. 0, or (b) to lplh and 2p2h configurations within second RPA approach in ref. |TT 



or (c) to three phonon configurations in ref. |12j , or (d) to all of them in ref. [13| . Somewhat 
aside from these two groups is the paper [[14]] where general properties of the DGDR are 
treated by means of a sum rule approach. 

In the present paper we investigate the contribution of non-natural parity 1 + two-phonon 
states to the total electromagnetic excitation cross sections. The [1~ <g> l~]i+ component of 
the DGDR was never considered in previous microscopic studies. They could, in principle, 
be responsible for the missing part of theoretical evaluations of the DGDR excitation cross 
sections. In phenomenological approaches describing the single GDR as one collective state, 
this component of the DGDR is forbidden by symmetry properties. Taking into account the 
Landau damping this collective state splits into a set of different 1^ states distributed over 
an energy interval, where i stands for the order number of each state. Again, the diagonal 
components [1^ ® ljf]i+ are forbidden by the same symmetry properties but nondiagonal 
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ones [17®1^7]i+, a priori, may be excited in two-step process bringing some "extra strength" 
in the DGDR region. Consequently, the role of these nondiagonal components depends on 
how strong is the Landau damping. 

To take into account the Landau damping effect we have performed an RPA calculation 
for J n = 1~ states in 208 Pb making use the quasiparticle-phonon model ||15|| . The model 
Hamiltonian includes an average field for protons and neutrons, treated by a Woods-Saxon 
potential, and a residual interaction in a separable form. For open shell nuclei it also has a 
term corresponding to monopole pairing. For our RPA calculation we have used the single 



particle spectrum and parameters of residual interaction from ref. [fLq| . This calculation 
provides us with the spectrum Ei of one-phonon 1~ states and reduced matrix elements 
< I7||-£71||<7.s. > of their electromagnetic excitation from the ground state. We produce 
two-phonon DGDR states with quantum numbers J 71 = + , 1 + and 2 + by coupling one- 
phonon RPA states with the wave function |17 > m , to each other. The index m stands for 
different magnetic substates. The wave function of the two-phonon states has the form: 

|[17 ® li~U=o+,2+ >m= -t= E {lmlm'\JM)\ir > m |17 > m ,, (0.1) 

v ^ m,m' 

for two-phonon states made of two identical phonons while for other DGDR states it is: 

|[17 ® l^=o+,i+, 2 + >m= E (lmlm'\JM)\l- > m (1,7 > m ,, (0.2) 

m,m' 

In the present calculation we do not include the interaction between DGDR states, of eqs. 
( P-l| - P.2|) , and we do not couple them to states with different than two number of phonons as 



was done in ref. |R|. Thus, our two-phonon states |[1~ <8> 17] >m have excitation energy 
equal to the sum of one-phonon energies Ei + Ey and are degenerated for different values of 
the total spin J n and its projection M. The reduced matrix element < [17 ® 17] \\E1 \ |17 > 
of electromagnetic excitation of two-phonon states, eqs. ( p.l| - p.2|) , from the one-phonon state 
1 17 > m is related, in the boson picture of nuclear excitation, to the excitation of |17 > m 
from the ground state as follows: 



< [17®17U||£1||17 >= W(1 + ^) (2J +1) < l«\\El\\g.s. > . (0.3) 



It should be noted that although for the two-phonon states, eq. ( |0.1| ), we have an extra factor 
\/2, the states of eq. ( p.2|) play a more important role in two-step excitations since they can 
be reached by two different possibilities: g.s. — > 17 — > [17 <8> 1.7/] and g.s. — > 1 4 7 — > [17 ® 1^7] . 
Making use of these nuclear structure ingredients we have performed calculations of DGDR 
excitation in relativistic heavy ion collisions in second-order perturbation theory and with 
a coupled-channels procedure. 
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The first effect we observed was that in second-order perturbation theory the amplitude 
for this process was identically zero in a semi-classical approach |T7[]. This can be under- 
stood by looking at figure |I|. The time-dependent field Ve\ carries angular momentum with 
projections m — 0, ±1. Thus, to reach the 1 + DGDR magnetic substates, many routes 
are possible. The lines represent transitions caused by the different projections of Vex'- (a) 
dashed lines are for m = 0, (b) dashed-dotted lines are for m = — 1, and (c) solid-lines 
are for m = +1. The relation VEi, m =o VEi,m=±i holds, so that not all routes yield the 
same excitation amplitude. Since the phases of the wavefunctions of each set of magnetic 
substates are equal, the difference between the transition amplitudes to a final M, can 
also arise from different values of the Clebsh-Gordan coefficients (lmlm'\lM). It is easy 
to see that, for any route to a final M, the second-order amplitude will be proportional 
to (001m|lm)(lmlm'|lM) VEi,m' Vei,™ + { m m ')- The two amplitudes carry opposite 
signs from the value of the Clebsh-Gordan coefficients. The identically zero result for the 
excitation amplitude of the 1 + DGDR state is therefore a consequence of 

^(00lTO|lm)(lmlm'|lM) = . (0.4) 

mm' 

We have also performed a coupled-channels calculation following the theory described 
in ref. [Q. As shown in ref. ||18|| , the coupling of the electric quadrupole (isovector and 



isoscalar) and the electric dipole states is very weak and can be neglected. We therefore 
include in our space only one-phonon 1~ and two-phonon [1 4 ~ £g> 1^]jt (J n = + , 1 + and 2 + ) 
states. We obtain the occupation amplitudes by solving the coupled-channels equations 



iha \ l " \ m =Y. < ^{ i- \ m \VEi(t)\V p, m > > exp ilEp-Et !- | 

l[l-<81-]J' m ~, l[l-®l-]J' m L V 



a/3, m > (0.5) 



where (3 = l^T or [1~ ® 1^]. The time dependent electric dipole field is that of a straight-line 
moving particle with charge Ze, and impact parameter b (we use eqs. (25-26) of ref. ||18|| ). 

Due to the large number of degenerate magnetic substates, to make our coupled-channels 
calculation feasible, we have chosen a limited set of GDR and DGDR states. We have taken 
six 1~ states which have the largest value of the reduced matrix element < l~||.E'l||<7.s. >. 
These six states exhaust 90.6 % of the classical EWSR, while all 1~ states up to 25 MeV in 
our RPA calculation exhaust 94.3 % of it. This value is somewhat smaller than the 122 % 



reported in ref. [19]. It is because the continuum in our RPA calculation was approximated 
by narrow quasibound states. From these six one-phonon 1~ states we construct two-phonon 
[l 4 r <S> 1jv]jt states, eqs. ( |0TT| - |0~2l ) , which also have the largest matrix element of excitation 
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< (g> lj7]jT||-El||lj r > for excitations starting from one-phonon states Q. The number 
of two-phonon states equals to twenty one for J 71 " = + and 2 + , and to fifteen for J w = 
1 + . The cross section for the DGDR excitation was obtained by a sum over the final 
magnetic substates of the square of the occupation amplitudes and, finally, by an integration 
over impact parameter. We have chosen the minimum impact parameter, b= 15.54 fm, 
corresponding to the parameterization of ref. ||21|| , appropriate for lead-lead collisions. 

The electromagnetic excitation cross sections for the reaction 208 Pb (640 MeV/nucleon) 
+ 208 Pb with excitation of all our basic 63 states is shown in figure |2|. The total cross 
sections for each multipolarity are presented in table |, together with the results of first-order 
(for one-phonon excitations) and second-order (for two-phonon excitations) perturbation 
theory. The coupled-channels calculation yields a non-zero cross section for the 1 + DGDR 
state due to other possible routes (higher-order), not included in second-order perturbation 
theory. One observes a considerable reduction of the DGDR cross sections, as compared to 
the predictions of the second-order perturbation theory. The GDR cross sections are also 
reduced in magnitude. However, the population of the 1 + DGDR states are not appreciable 
and cannot be the source of the missing excitation cross section needed to explain the 
experiments. In general, the coupled-channels calculation practically does not change the 
relative contribution of different one-phonon 1~ and two-phonon states [1^ <8> to the 

total cross section with given J 71 = 1~, + and 2 + . But since the 1 + component of DGDR, 
with its zero value of excitation cross section in the second-order perturbation theory, has 
a special place among the two other components, the main effect of coupled-channels is to 
redistribute the total cross section between the J n = + , 2 + and J n = 1 + components. 

The calculated cross section in coupled-channels for both GDR and DGDR are somewhat 
smaller than reported in experimental findings H . This is not surprising since as mentioned 
above our chosen six 1~ states exhaust only 90.6 % of EWSR while the photo-neutron data 



l lt was demonstrated in ref. JT^,f20| that direct excitation of two-phonon configurations from the 
ground state is very weak. It allows us to exclude in our calculation matrix elements of the form 
< [1^ (g) l i 7] 2 +(i+)||i?2(Ml)||g.s. > which correspond to direct transitions and produce higher order 
effects in comparison with accounted ones. These matrix elements give rise to DGDR excitation 
in first order perturbation theory. Thus, to prove our approximation we have calculated such cross 
sections and got total values equal to 0.11 mb and jO.Ol mb for the twenty one 2 + and the fifteen 
1 + basic two-phonon states, respectively. These values have to be compared to 244.9 mb for the 
total DGDR cross section in the second order perturbation theory. 
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T9fl indicate that this value equals to 122 %. Due to this underestimate of exhaust of the 



EWSR the cross section for DGDR excitation reduces more strongly than the one for the 
single GDR. This is because the GDR cross section is roughly proportional to the total 
B(E1) value while for DGDR it is proportional to the square of it. If we apply a primitive 
scaling to obtain the experimental value 122 % of EWSR the ratio R = ct(dgdk)/ 'v(gdr), 
the last line of our table |, changes into 0.096 and 0.101 for the coupled-channels calculation 
and for the perturbation theory, respectively. The experimental findings || yield the value 
Rexp = 0.116±0.014. The reported || disagreement R exp / R ca ic =1.33±0.16 is the result of 
a comparison with R ca ic obtained within a folding model, assuming 122 % of the EWSR. 
We get a somewhat larger value of R ca ic (taking into account our scaling procedure) because 
the B(E1) strength distribution over our six 1~ states is not symmetrical with respect to 
the centroid energy, Egdr'- the lower part is enhanced. A weak energy dependence in the 
excitation amplitude (see also ref. |§ for a discussion of this problem), which is also squared 
for the DGDR, enhances the DGDR cross section for a non-symmetrical distribution with 
respect to the symmetrical one, or when the GDR is treated as a single state. The effect of 
the energy dependence is demonstrated for a single GDR in the top part of figure |2| where 
the excitation cross sections are compared to the B(E1) strength distribution. It produces a 
shift to lower energies of the centroid of the GDR and the DGDR cross sections with respect 
to the centroid of the B(E1) and the B(El)xB(El) strength distribution, respectively. In 
our calculation this shift equals to 0.26 MeV for the GDR and to 0.33, 0.28 MeV for the 
DGDR within coupled-channels and perturbation theory, respectively. 

Of course, this scaling procedure has no deep physical meaning but we have included 
this discussion to indicate that the disagreement between experiment and theory for the 
DGDR excitation cross sections in 208 Pb reached the stage when theoretical calculations 
have to provide a very precise description of both the GDR and the DGDR to draw up final 
conclusions. Work in this direction is in progress now and will be reported soon. 

In summary, we have studied the excitation cross sections of the 1 + component of the two- 
phonon giant dipole resonance in 208 Pb (together with + and 2 + components) in relativistic 
heavy ion collisions. In second order perturbation theory this cross section equals to zero 
due to angular momentum properties. In coupled-channels calculations this component 
contributes to the total cross section. But it is appreciably quenched with respect to other 
components of DGDR. We indicate that a precise microscopic calculation of the giant dipole 
resonance is required to answer the question whether the problem of enhanced experimental 
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cross section for DGDR excitation in relativistic heavy ion collisions stands for or not. 
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TABLES 

TABLE I. Cross section (in mb) for the excitation of the GDR and the three components with 
= 0+, 2+, 1+ of the DGDR in 208 Pb (640 MeV/A) + 208 Pb collisions. Calculations are performed 

within coupled-channels (CC) and within the first (PT-1) and second (PT-2) order perturbation 

theory, respectively. 





CC 


PT-1 


PT-2 


GDR 


2830. 


3275. 


0. 


DGDR + 


33.0 


0. 




43.1 


DGDR 2 + 


163.0 


0.11 




201.8 


DGDR l+ 


6.3 


jO.Ol 




0. 


DGDR/GDR 


0.071 


0.075 
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FIGURES 

FIG. 1. The possible paths to the excitation of a given magnetic substate of the 1 + component 
of the DGDR are displayed. The transitions caused by the different projections of the operator Vei 
are shown by: (a) dashed lines for m = 0, (b) dashed-dotted lines for m = — 1, and (c) solid-lines 
are m = +1. 

FIG. 2. The electromagnetic excitation cross sections for the reaction 208 Pb (640 MeV/nucleon) 
+ 208 Pb calculated in coupled-channels. It is shown the excitation of the GDR (top) and the three 
components J 71 arbitrary units) over 1 _ states is shown by dashed lines. For visuality it is shifted 
up by 100 keV. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9606004v 1 



This figure "fig2-l.png" is available in "png" format from: 
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